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Abstract 

In this paper, we investigate the asymptotic spectrum of complex or real Deformed Wigner 
matrices {Mn)n defined by Afjv = Wn/VT^ + An where Wn is a N x N Hermitian (resp. 
symmetric) Wigner matrix whose entries have a symmetric law satisfying a Poincare inequality. 
The matrix An is Hermitian (resp. symmetric) and deterministic with all but finitely many 
eigenvalues equal to zero. We first show that, as soon as the first largest or last smallest eigenvalues 
of An axe sufficiently far from zero, the corresponding eigenvalues of Mn almost surely exit the 
limiting semicircle compact support as the size becomes large. The corresponding limits are 
universal in the sense that they only involve the variance of the entries of Wn- On the other hand, 
when An is diagonal with a sole simple non-null (fixed) eigenvalue large enough, we prove that the 
fluctuations of the largest eigenvalue are not universal and vary with the particular distribution 
of the entries of Wjv. 

1 Introduction 

This paper lies in the Uneage of recent works studying the influence of some perturbations on the 
asymptotic spectrum of classical random matrix models. Such questions come from Statistics (cf. 
[Jo]) and appeared in the framework of empirical covariance matrices, also called non- white Wishart 
matrices or spiked population models, considered by J. Baik, G. Ben Arous and S. Peche [Bk-B-P] 
and by J. Baik and J. Silverstein [Bk-Sl]. The work [Bk-B-P] deals with random sample covariance 
matrices {Sn)n defined by 

Sn = j^Y^Y^ (1.1) 

where Y/v is a p x iV complex matrix whose sample column vectors are i.i.d, centered, Gaussian and of 
covariance matrix a deterministic Hermitian matrix Ep having all but finitely many eigenvalues equal 
to one. Besides, the size of the samples N and the size of the population p = pM are assumed of the 
same order (as ^ oo). The authors of [Bk-B-P] first noticed that, as in the classical case (known 
as the Wishart model) where Sp = Ip is the identity matrix, the global limiting behavior of the spec- 
trum of Sn is not affected by the matrix Sp. Thus, the limiting spectral measure is the well-known 
Marchenko-Pastur law. On the other hand, they pointed out a phase transition phenomenon for the 
fluctuations of the largest eigenvalue according to the value of the largest eigenvalue (s) of Sp. The 
approach of [Bk-B-P] does not extend to the real Gaussian setting and the whole analog of their result 
is still an open question. Nevertheless, D. Paul was able to establish in [P] the Gaussian fluctuations 
of the largest eigenvalue of the real Gaussian matrix Sm when the largest eigenvalue of Ep is simple 
and sufficiently larger than one. More recently, J. Baik and J. Silverstein investigated in [Bk-Sl] the 
almost sure limiting behavior of the extremal eigenvalues of complex or real non necessarily Gaussian 
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matrices. Under assumptions on the first four moments of the entries of Yn, they showed in particu- 
lar that when exactly k eigenvalues of Sp are far from one, the k first eigenvalues of Sn are almost 
surely OTitside the limiting Marchenko-Pastur support. Fluctuations of the eigenvalues that jump are 
universal and have been recently found by Z. Bai and J. F. Yao in [B-Ya2] (we refer the reader to 
[B-Ya2] for the precise restrictions made on the definition of the covariance matrix Sp). Note that 
the problem of the fluctuations in the very general setting of [Bk-Sl] is still open. 

Our purpose here is to investigate the asymptotic behavior of the first extremal eigenvalues of 
some complex or real Deformed Wigner matrices. These models can be seen as the additive analogue 
of the spiked population models and are defined by a sequence (MAr)jv given by 

Mn = -^Wn + An := Xn + An (1.2) 

where Wn is a Wigner matrix such that the common distribution of its entries satisfied some tech- 
nical conditions (given in (i) below) and is a deterministic matrix of finite rank. We establish 
the analog of the main result of [Bk-Sl] namely that, once An has exactly k (fixed) eigenvalues far 
enough from zero, the k first eigenvalues of Mjv jump almost surely outside the limiting semicircle 
support. This result is universal (as the one of [Bk-Sl]) since the corresponding almost sure limits 
only involve the variance of the entries of Wn- On the other hand, at the level of the fiuctuations, 
we exhibit a striking phenomenon in the particular case where An is diagonal with a sole simple 
non-null eigenvalue large enough. Indeed, we find that in this case, the fluctuations of the largest 
eigenvalue of Mn are not universal and strongly depend on the particular law of the entries of Wn- 
More precisely, we prove that the limiting distribution of the (properly rescaled) largest eigenvalue of 
Mn is the convolution of the distribution of the entries of Wn with a Gaussian law. In particular, if 
the entries of Wn are not Gaussian, the fluctuations of the largest eigenvalue of Mn are not Gaussian. 

In the following section, we first give the precise definition of the Deformed Wigner matrices (1.2) 
considered in this paper and we recall the known results on their asymptotic spectrum. Then, we 
present our results and sketch their proof. We also outline the organization of the paper. 



2 Model and results 

Throughout this paper, we consider complex or real Deformed Wigner matrices {Mn)n of the form 
(1.2) where the matrices Wn and An are defined as follows: 

(i) Wn is a. N X N Wigner Hermitian (resp. symmetric) matrix such that the N"^ random variables 
{Wn)u, V^^e{{WN)ij)i<j, ^/2Qm{{WN)ij)i<j (resp. the random variables -^{WN)ii, 
{WN)ij, i < j) are independent identically distributed with a symmetric distribution ^ of vari- 
ance (T^ and satisfying a Poincare inequality (see Section 3). 

(ii) is a deterministic Hermitian (resp. symmetric) matrix of fixed finite rank r and built from 
a family of J fixed real numbers 9i > ■ ■ ■ > 9 j independent of N with some jo such that 
^3o — 0- We assume that the non-null eigenvalues 9j oi An are of fixed multiplicity kj (with 
T^j^j^ kj = r) i.e. Ajv is similar to the diagonal matrix 

Dn = diag(6i, . . . ,6i,. . . , 6j i, . . . , 6j ,0, Oj +i, . . . , Oj j^i, . . . ,6j,. . . , 6j). (2.1) 

^ V ' ^ V ' ' V ' ^ ^ ' ^ V ' 

Before going into details of the results, we want to point out that the condition made on (namely 
that /i satisfies a Poincare inequality) is just a technical condition: we conjecture that our results still 
hold under weaker assumptions (see Remark 2.1 below). Nevertheless, a lot of measures satisfy a 
Poincare inequality (we refer the reader to [B-G] for a characterization of such measures on M, see 
also [A. et al]). For instance, consider ^{dx) = exp(— |x|")rfa; with a>l. 

Futhermore, note that this condition implies that fi has moments of any order (cf. Corollary 3.2 and 
Proposition 1.10 in [L]). 



Let us now introduce some notations. When the entries of Wn are further assumed to be Gaussian 
that is, in the complex (resp. real) setting when Wn is of the so-called GUE (resp. GOE), we will 
write instead of Wn- Then X§ := -^W^ will be said of the GU(0)E(iV, ^) and we will let 
M§ = X§ + An he the corresponding Deformed GU(0)E model. 

In the following, given an arbitrary Hermitian matrix B of order N, we will denote by Ai (B) > ■ ■ ■ > 
Xn{B) its A'' ordered eigenvalues and by hb = jf^i=i^>^i(B) empirical measure. For notational 
convenience, we will also define Ao(-B) = +00 and Ajv+i(-B) = —00. 

The Deformed Wigncr model is built in such a way that the Wigner Theorem is still satisfied. 
Thus, as in the classical Wigner model {An = 0), the spectral measure {hmn) converges a.s. to the 
semicircle law Hsc whose density is given by 



dx ^""^ = ^V4a2-^^l[-2.,2.](a=). (2.2) 

This result follows from Lemma 2.2 of [B]. Note that it only relies on the two first moment assump- 
tions on the entries of Wn and the fact that the An^s are of finite rank. 

On the other hand, the asymptotic behavior of the extremal eigenvalues may be affected by the 
perturbation An ■ Recently, S. Pcchc studied in [Pe] the Deformed GUE under a finite rank perturba- 
tion An defined by (ii). Following the method of [Bk-B-P], she highlighted the effects of the non-null 
eigenvalues oi An at the level of the fluctuations of the largest eigenvalue of M^. To explain this in 
more detail, let us recall that when An = 0, it was established in [T-W] that as A?' — > 00, 

a-^N^"'(\i{X%)-2a) ^ F2 (2.3) 

where F2 is the well-known GUE Tracy- Widom distribution (see [T-W] for the precise definition). 
Dealing with the Deformed GUE M^, it appears that this result is modified as soon as the first 
largest eigenvalue(s) of An are quite far from zero. In the particular case of a rank one perturbation 
An having a fixed non-null eigenvalue ^ > 0, [Pe] proved that the fluctuations of the largest eigenvalue 
of are still given by (2.3) when 6 is small enough and precisely when 6 < a. The limiting law is 
changed when 6 = a. As soon as ^ > cr, [Pe] established that the largest eigenvalue Ai(M^) fluctuates 
around 
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Pe = + ^ (2.4) 

(which is > 2fT since 6 > a) as 

\/iv(Ai(M^) - pe) ^(0, al) (2.5) 

where 



Similar results are conjectured for the Deformed GOE but S. Peche emphasized that her approach 
fails in the real framework. Indeed, it is based on the explicit Predholm determinantal representation 
for the distribution of the largest eigenvalue(s) that is specific to the complex setting. Nevertheless, 
M. Maida [M] obtained a large deviation principle for the largest eigenvalue of the Deformed GOE 
M§ under a rank one deformation An; from this result she could deduce the almost sure limit 
with respect to the non-null eigenvalue oi An- Thus, under a rank one perturbation An such that 
Dn = diag(6', 0, • • • ,0) where 6* > 0, [M] showed that 

Ai(M^)^p,, iie>a (2.7) 

and 

Ai(M^) ^ 2(7, if ^ < (T. (2.8) 



Note that the approach of [M] extends with minor modifications to tlic Deformed GUE. Following the 
investigations of [Bk-Sl] in the context of general spiked population models, one can conjecture that 
such a phenomenon holds in a more general and non necessarily Gaussian setting. The first result of 
our paper, namely the following Theorem 2.1, is related to this question. Before being more explicit, 
let us recall that when = 0, the whole spectrum of the rescaled complex or real Wigner matrix 
= "^^w belongs almost surely to the semicircle support [— 20-, 2cr] as N goes to infinity and 
that (cf. [B-Yi] or Theorem 2.12 in [B]) 

Xi{Xn) 2a and Xn{Xn) ^ -2a. (2.9) 

Note that this last result holds true in a more general setting than the one considered here (see [B-Yi] 
for details) and in particular only requires the finiteness of the fourth moment of the law fi. Moreover, 
one can readily extend the previous limits to the first extremal eigenvalues of Xj^ i.e. 

for any fixed A; > 1, Xk{X]\[) 2a and XN-k{XN) — — * —2a. (2.10) 

Here, we prove that, under the assumptions (i)-(ii), (2.10) fails when some of the 9j's are sufficiently 
far from zero: as soon as some of the first largest (resp. the last smallest) non- null eigenvalues 6j 's of 
Apf are taken strictly larger then a (resp. strictly smaller than —a), the same part of the spectrum of 
Mjv almost surely exits the semicircle support [—2a, 2a] as AT — > oo and the new limits are the pe^ 's 
defined by 
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/'%=^i + ^- (2-11) 

Observe that pe- is > 2a (resp. < —2a) when 0j > a (resp. < —a) (and p0. = ±2a if 9j = ±cr). 
Here is the precise formulation of our result. For definiteness, we set fci H + kj-\ := if j = 1. 

Theorem 2.1. Let J+a (resp. he the number of j's such that 0j > a (resp. dj < —a). 

(a) VI < j < J+a, VI < i < fcj, Xki+...+kj-1+iiMN) — > PBj a.s., 

(b) Xk^+...+kj^^+i{MN) — > 2(7 a.s., 

(c) Xk^+...+kj-j_^{MN) — ^-2a a.s., 

(d) Vj > J- J_<^ + 1, VI < i < fcj, Afei+...+fe^._i+i(Mjv) — > P0J a.s. 



Remark 2.1. Let us notice that, following [Bk-Sl], one can expect that this theorem holds true in a 
more general setting than the one considered here, namely only requires four first moment conditions 
on the law p of the Wigner entries. As we will explain in the following, the assumption (i) that 
p satisfies a Poincare inequality is actually fundamental in our reasoning since we will need several 
variance estimates. 

This theorem will be proved in Section 4. The second part of this work is devoted to the study of 
the particular rank one diagonal deformation An = diag(0, 0, • • • ,0) such that 9 > a. We investigate 
the fluctuations of the largest eigenvalue of any real or complex Deformed model Mat satisfying (i) 
around its limit pg (given by the previous theorem). Wc obtain the following result. 

Theorem 2.2. Let A^ = diag(0, 0, • • • ,0) and assume that 9 > a. Define 

^^=4(^^) + 2^ (2-12) 

where 1 = 4 (resp. t = 2) when Wn is real (resp. complex) and m4 := / x'^dfj.{x). Then 

2 

y/N(^Xi{MN)-pe) {1 ^ ^){p*Af{0,vg)]. (2.13) 



Note that when = 2>a^ as in the Gaussian case, the variance of the hmiting distribution in (2.13) 
is equal to (resp. 2cr^) in the complex (resp. real) setting (with ae given by (2.6)). 

Remark 2.2. Since /x is symmetric, it readily follows from Theorem 2.2 that when An = diag(^, 0, • • • , 0) 
and 9 < —a, the smallest eigenvalue of fluctuates as 



In particular, one derives the analog of (2.5) for the Deformed GOE that is 

Theorem 2.3. Let An he an arbitrary deterministic symmetric matrix of rank one having a non-null 
eigenvalue 9 such that 9 > a. Then the largest eigenvalue of the Deformed GOE fluctuates as 



Obviously, thanks to the orthogonal invariance of the GOE, this result is a direct consequence of 
Theorem 2.2. 

It is worth noticing that, according to the Cramer-Levy Theorem (cf. [F], Theorem 1 p. 525), 
the limiting distribution (2.13) is not Gaussian if fi is not Gaussian. Thus, (2.13) depends on the 
particular law /x of the entries of the Wigner matrix Wn which implies the non-universality of the 
fluctuations of the largest eigenvalue of rank one diagonal deformation of symmetric or Hermitian 
Wigner matrices (as conjectured in Remark 1.7 of [Fc-Pc]). 

The latter also shows that in the non-Gaussian setting, the fluctuations of the largest eigenvalue 
depend, not only on the spectrum of the deformation An, but also on the particular definition of 
the matrix An- Indeed, in collaboration with S. Pcchc, the third author of the present article has 
recently stated in [Fe-Pe] the universality of the fluctuations of some Deformed Wigner models under 
a full deformation An deflned by (^Ar)ij = for all 1 < < A'' (see also [Fu-K]). Before giving 
some details on this work, we have to precise that [Fc-Pe] considered Deformed models such that the 
entries of the Wigner matrix Wn have sub-Gaussian moments. Nevertheless, thanks to the analysis 
made in [R], one can observe that the assumptions of [Fe-Pe] can be reduced and that it is for example 
sufficient to assume that the Wij's have moment of order 19 (the precise condition of [R] is given 
by (2.15) below). Thus, the conclusions of [Fe-Pe] apply to the setting considered in our paper. The 
main result of [Fc-Pc] establishes the universality of the fluctuations of the largest eigenvalue of the 
complex Deformed model Mn associated to a full deformation An and for any value of the parameter 
9. In particular, when 9 > a, it is proved therein the universality of the Gaussian fluctuations (2.5). 
Notice that the approach of [Fe-Pe] is completely different from the one developed below in Section 
5 to derive Theorem 2.2. It is mainly based on a combinatorial method inspired by the work [So] 
(which handles the non-Deformed Wigner model) and the known fluctuations for the Deformed GUE 
(given by [Pe]). The combinatorial arguments of [Fe-Pe] also work (with minor modifications) in the 
real framework and yields the universality of the fluctuations if 9 < a. In the case where 9 > a which 
is of particular interest here, the analysis made in [Fe-Pe] reduced the universality problem in the real 
setting to the knowledge of the particular Deformed GOE model which was unknown up to now (note 
that this remark is also valid in the case where 9 = a). Thus, thanks to our previous Theorem 2.3 
and the analysis of [Fe-Pe] and [R], we are now in position to claim the following universality. 

Theorem 2.4. Let An he a full perturbation given by {An)^ = jf for all {i,j)- Assume that 9 > a. 
Let Wn be an arbitrary real Wigner matrix with the underlying measure fi being symmetric with a 
variance and such that there exists some p> 18 satisfying 



Then the largest eigenvalue of the Deformed model Mn has the Gaussian fluctuations (2.14)- 

Remeirk 2.3. To be complete, let us notice that the previous result still holds when we allow the 
distribution v of the diagonal entries of Wn being different from n provided that v is symmetric and 
satisfies (2.15). 





(2.14) 



(2.15) 



The fundamental tool of this paper is the Stieltjcs transform. For z € C\M, we denote the resolvent 
of the matrix M jv by 

G,vW = (2/iv-Mw)-i 

and the Stieltjes transform of the expectation of the empirical measure of the eigenvalues of Mjv by 

5jv(-2) =E(trjv(Gjv(^))) 
where trjv is the normalized trace. We also denote by 

the Stieltjes transform^ of a random variable s with semi-circular distribution 

Theorem 2.1 is the analog of the main statement of [Bk-Sl] established in the context of general 
spiked population models. The conclusion of [Bk-Sl] requires numerous results obtained previously 
by J. Silverstein and co-authors in [Ch-S], [B-Sl] and [B-S2] (a summary of all this literature can be 
found in [B] pp. 671 675). From very clever and tedious manipulations of some Stieltjes transforms 
and the use of the matricial representation (1.1), these works highligh a very close link between the 
spectra of the Wishart matrices and the covariance matrix (for quite general covariance matrix which 
include the spiked population model). Our approach mimics the one of [Bk-Sl]. Thus, using the 
fact that the Deformed Wigner model is the additive analog of the spiked population model, several 
arguments can be quite easily adapted here (this point has been explained in Chapter 4 of the PhD 
Thesis [Fe]). Actually, the main point in the proof consists in establishing that for any £ > 0, almost 
surely, 

Spect(M;v)ci^|(ei,--- ,ej) (2.16) 
for all N large, where we have defined 

and 

,0j) := \^Pej;Pej-i;--- ;Pej-j_„+i}u[-2(7;2a]u|pej+^;pej+^_j;--- ;pe^Y 

This point is the analog of the main result of [B-Sl]. The analysis of [B-Sl] is based on technical and 
numerous considerations of Stieltjes transforms strongly related to the Wishart context and that can 
not be directly transposed here. Thus, our approach to prove such an inclusion of the spectrum of 
Mtv is very different from the one of [B-Sl]. Indeed, we use the methods developed by U. Haagerup 
and S. Thorbj0rnsen in [H-T], by H. Schultz [S] and by the two first authors of the present article 
[C-D] . The key point of this approach is to obtain a precise estimation at any point z G C\M of the 
following type 

g,{z)-gNiz) + j^L,{z) = 0{^), (2.17) 

where L^- is the Stieltjcs transform of a distribution Ao- with compact support in Ku{6i, ■ ■ ■ , 9j). In- 
deed such an estimation allows us through the inverse Stieltjes transform and some variance estimates 
to deduce that trAr lcxe(ei,. . .0j)(Mjv) = 0{1/N3) a.s.. Thus the number of eigenvalues of Mjv in 
''K^{6i, ■ ■ ■ , 9j) is almost surely a 0(1/A^^) and since for each A'' this number has to be an integer, 
we deduce that it is actually equal to zero as N goes to infinity. 

Dealing with the particular diagonal perturbation = diag(0, 0, . . . , 0), we obtain the fluctua- 
tions of the largest eigenvalue Xi^Mn) by an approach close to the one of [P] and the ideas of [B-B-P]. 
The reasoning relies on the writing of Ai(Mjv) in terms of the resolvent of a non-Deformed Wigner 
matrix. 

The paper is organized as follows. In Section 3, we introduce preliminary lemmas which will be of 
basic use later on. Section 4 is devoted to the proof of Theorem 2.1. We first establish an equation 
(called "Master equation") satisfied by gN up to some correction of order (see Section 4.1). Then 



^Note that in some papers for which we make reference, the Stieltjes transform is defined with the opposite sign. 



wc explain how this master equation gives rise to an estimation of type (2.17) and thus to the inelusion 
(2.16) of the spectrum of Mm in K^{di, ■ ■ ■ , 6j) (sec Sections 4.2 and 4.3). In Section 4.4, we use this 
inclusion to relate the asymptotic spectra oi Am and Mjv and then deduce Theorem 2.1. The last 
section states Theorem 2.2. 

Acknowledgments. The authors are very grateful to Jack Silverstein and Jinho Baik for providing 
them their proof of Theorem 5.2 (that is a fundamental argument in the proof of Theorem 2.2) which 
is presented in the Appendix of the present article. 

3 Basic lemmas 

We assume that the distribution fi of the entries of the Wigner matrix Wn satisfies a Poincare 
inequality: there exists a positive constant C such that for any function / : M ^ C such that / 

and /' are in L'^{iJ.), 

v(/) <cj \ffdn, 

withV(/)=E(|/-E(/)|2). 

For any matrix M, define ||M||2 = (Tr(M*M))i/2 the Hilbert-Schmidt norm. Let ^ : (MAr(C)sa) ^ 
(resp. : (Mjv(C)s) — > K ' ' ) be the canonical isomorphism which maps an Hermitian (resp. 
symmetric) matrix to the real parts and the imaginary parts of its entries (resp. to the entries) 

iM)i,,i<j. 

Lemma 3.1. Let Mn be the complex (resp. real) Wigner Deformed matrix introduced in Section 

2. For any function f : {resp. ]S.~^ ~) C such that f and the gradient V(/) are both 
polynomially bounded, 

V[f o *(M;v)] < ^E{|| V [/ o ^r(Mjv)] 11^}. (3.1) 
Proof: According to Lemma 3.2 in [C-D], 

V[/ o *(X;v)] < ^E{|| V [/ o *(X;v)] 11^}. (3.2) 

Note that even if the result in [C-D] is stated in the Hermitian case, the proof is valid and the result 

still holds in the symmetric case. Now (3.1) follows putting g{xij\i < j) := f{xij + (^jv)ij;i < j) in 
(3.2) and noticing that the {AN)ij are uniformly bounded in i,j,N. □ 

This lemma will be useful to estimate many variances. 

Now, we recall some useful properties of the resolvent (see [K-K-P], [C-D]). For any Hermitian matrix 
M we denote its spectrum by Spect(M). 

Lemma 3.2. For a N x N Hermitian or symmetric matrix M, for any z G C\Spect(M), we denote 
by G{z) := (z/jv - M)~^ the resolvent of M. 

Let z e C\M, 

(i) ||G(-2;)|| < |9m(z)|~^ where ||.|| denotes the operator norm. 

(ii) \G{z)ij\ < \'^m{z)\-'^ for all i,j = 1,...N. 
(Hi) For p > 2, 

(iv) The derivative with respect to M of the resolvent G{z) satisfies: 



G'j^{z).B = G{z)BG{z) for any matrix B. 



(v) Let z e C such that \z\ > ||M||; we have 



\\G{z)\\ < 



M\ 



Proof: We just mention that (v) comes readily noticing that the eigenvalues of the normal matrix 
G{z) are the = 1, . . . ,N.n 

We will also need the following estimations on the Stieltjes transform g^r of the semi-circular distribu- 
tion /Use- 

Lemma 3.3. is analytic on <C\[—2a,2a] and 
• Vz G {z G C : Qmz 0}, 

a^gliz) - zg^iz) + 1 = 0. (3.3) 
\g^{z)\ < \Qmz\-\ (3.4) 
\9cr{z)~'^\ < \z\ + a'^\Qmz\-^. (3.5) 



For a > 0,6 € M, 
\/z€{zeC: \z\> 2a}, 



{z-t) 

1 



aga{z) - z + 9 



< \'^m{z)\-\ (3.7) 



M^)|-^<N + ^. (3.10) 

Proof: For the equation (3.3), we refer the reader to Section 3.1 of [B]. (3.7) is a consequence of 
?sm{gcr{z))?sm{z) < 0. Other inequalities derive from (3.3) and the definition of ga- □ 

4 Almost sure convergence of the first extremal eigenvalues 

Sections 4.1, 4.2 and 4.3 below describe the different steps of the proof of the inclusion (2.16). We 

choose to develop the case of the complex Deformed Wigner model and just to point out some differ- 
ences with the real model case (at the end of Section 4.3) since the approach would be basically the 
same. In these sections, we will often refer the reader to the paper [C-D] where the authors deal with 
several independent non Deformed Wigner matrices. The reader needs to fix r = 1, m = 1, oq = 0, 
ai = a and to change the notations X = z, Gn = gN, G = jfo- in [C-D] in order to use the different 
proofs we refer to in the present framework. We shall denote by Pk any polynomial of degree k with 
positive coefficients and by C, K any constants; Pk, C, K can depend on the fixed eigenvalues of 
and may vary from line to line. We also adopt the following convention to simplify the writing: we 
sometimes state in the proofs below that a quantity An{z), z G C\M is 0{^), p= 1,2. This means 
precisely that 



for some k and some I and we give the precise maj oration in the statements of the theorems or 
propositions. 

Section 4.4 explains how to deduce Theorem 2.1 from the inclusion (2.16). 



4.1 The master equation 
4.1.1 A first master inequality 

In order to obtain a master equation for gN{z), we first consider the Gaussian case, i.e. Xj^^ — ^® 
distributed as tlie GUE(iV, ^) distribution.^ 

Let us recall the integration by part formula for the Gaussian distribution. 

Lemma 4.1. Let ^ be a complex valued function on (M;v(C)sa) and Xn ~ GUE(N, ^ ). Then, 

N 

E[ct>'{XN).H] = —E[<P{Xn) Tr{XNH)] (4.1) 

for any Hermitian matrix H, or by linearity for H = Ejk, 1 < j,k < N where Ejk, 1 < j,k < N is 
the canonical basis of the complex space of N x N matrices. 

We apply the above lemma to the fimction ^{Xjs[) = {GN{z))ij = {{zl^ — Xjq — Ais[)^^)ij, 
z e C\M, 1 < i,j < N. In order to simplify the notation, we write {GN{z))ij — Gij. We obtain, for 
H = Eij : 

E{{GHG)ij) = ^E{GijTT{XNH)] 
N 

E{GuGjj) = -^E{Gij{XN)ji] 
Now, we consider the normalized sum of the previous identities to obtain: 

E((tr;vG)') = 4lE(trjv(GX;v)). 

(7 

Then, since 

GXn = {z-Xn- An^^Xn +An- zIn - An + zIn) = -In - GAn + zG, 
we obtain the following master equation: 

E((trAr Gf) + ^(-^E(trAr G) + 1 + E(trAr GAn)) = 0. 

Now, it is well known (see [C-D], [H-T] and Lemma 3.1) that: 

Var(trjv(G)) < ' ' 

thus, we obtain: 

Proposition 4.1. The Stieltjes transform satisfies the following inequality: 

\a^gl{z) - zgN{z) + 1 + ^E{T^{Gn{z)An))\ < C^^^ (4.2) 

Note that since An is of finite rank, E(Tr(Giv(^)^Jv)) < C where C is a constant independent of N 
(depending on the eigenvalues of An and z). 

We now explain how to obtain the corresponding equation (4.2) in the Wigner case. Since the 
computations are the same as in [C-D]^ and [K-K-P]^, we just give some hints of the proof. 

Step 1: The integration by part formula for the Gaussian distribution is replaced by the following 
tool: 



^Throughout this section, we will drop the subscript G in the interest of clcirity. 
^This paper treats the case of several independent non Deformed Wigner matrices. 
''The authors considered a non Deformed Wigner matrix in the symmetric real setting. 



Lemma 4.2. Let £, be a real-valued random variable such that E(|^|^'+^) < oo. Let 4> be a function 
from R to C such that the first p+1 derivatives are continuous and bounded. Then, 



(4.3) 



a=0 



where Ha are the cumulants o/^, |e| < Csup^ (<)|E(|^|P"'"^), C depends on p only. 



We apply this lemma with the function (/>(^) given, as before, by 4>{^) = Gij and ^ is now one of 
the variable 3?e((Xjv)fei)) ^"^((-'^Jv)fc;)- Note that, since the above random variables are symmetric, 
only the odd derivatives in (4.3) give a non null term. Moreover, as we arc concerned by estimation of 
order of (?jv, we only need to consider (4.3) up to the third derivative (see [C-D]). The computation 
of the first derivative will provide the same term as in the Gaussian case. 

Step 2: Study of the third derivative. 

We refer to [C-D] or [K-K-P] for a detailed study of the third derivative. Using some bounds on 
GjV) see Lemma 3.2, we can prove that the only term arising from the third derivative in the master 
equation, giving a contribution of order is: 



1e 

N 



fc=i 



In conclusion, the first master equation in the Wigner case reads as follows: 
Theorem 4.1. For z G C\K, gwiz) satisfies 



1 



1 K4 



a'gNizY - ZQNiz) + 1 + -mTr{GNiz)Ar,)] + -^E 



N 



where K4 is the fourth cumulant of the distribution /U. 



N 2 



< 



N 



l^(G^(z))^ 

fc=i 

P6(|3m(2)|-i) 



N2 



(4.4) 



4.1.2 Estimation of j^jv — ga] 

Since 

|E[Tr(G;v(^)^jv)]| V |E 
Theorem 4.1 implies that for any z g 



' 1 

-J2iGr,iz))l 



k=l 



<P4\Qm{z)\-'), 



\a^gN{zf - zgN{z) + 1| < 



P6(|3m(z)|-i) 



TV 



(4.5) 



To estimate \gN — ffo-| from the equation (3.3) satisfied by the Stieltjes transform g„ on the one hand 
and from the equation (4.5) on the other hand, we follow the method initiated in [H-T] and [S]. We 
don't develop it here since it follows exactly the lines of Section 3.4 in [C-D] but wc briefly recall the 
main arguments and results which will be useful later on. We define the open connected set 



0'^ = {z& C, ^m{z) > 0, :^&^(a2|SJm(^)|-i + |^|) < ^ 



4|3mz| 



One can prove that for any z in Oj^, 
• gN^z) ^ and 



\gN{z) 



<2{a^\'^m{z)\-^ + \z\) 



(4.6) 



An{z) := a^gN^z) + is such that 



and we have 



\Am{z) -z\< M^^i^2{a'\Qmz\-' + \z\) (4.7) 



9m(A^(z))>5^>0. (4.8) 



• Writing the equation (3.3) at the point An{z), we easily get that 

gN{z) = ga{AN{z)) (4.9) 

on the non empty open subset 0"j^ = {z G O'^f, '^mz > \/2a} and then on O'j^ by the principle 
of uniqueness of continuation. 

This allows us to get an estimation of \gN{z) — ga{z)\ on Oj^ and then to deduce: 

Proposition 4.2. For any z gC such that Qm{z) > 0, 

\gN{z) - 9Az)\ < {\z\ + i^)^^^^^. (4.10) 
4.1.3 Study of the additional term E[Tr(AjvGjv(^))] 

Prom now on and until the end of Section 4.1, we denote by 71, . . . , 7^ the non- null eigenvalues of Ajv 

{li = for some j =/= jo) in order to simplify the writing. Lot Un := U he a unitary matrix such that 
An = U*AU where A is the diagonal matrix with entries An = 7^, i < r; An =0, i > r. We set 

hN{z)=E[Ti{ANGN{z))] (4.11) 

r N 
fe=l i,j=l 

Our aim is to express hN{z) in terms of the Stieltjes transform gniz) for A'' large, using the integration 
by part formula. Note that since we want an estimation of order 0{^) in the master inequality (4.4), 
we only need an estimation of h]\[{z) of order 0{j^). As in the previous subsection, we first write the 
equation in the Gaussian case and then study the additional term (third derivative) in the Wigner 
case. 

a) Gaussian case 
Apply the formula (4.1) to $(Xjv) = Gji and H = En to get 

N 

E[GjiGii] = —^E[Gji{XN)ii] 
a 

and 

N 



1^E[G,,G„] = -^E[(GXa,),,]. 



N 

Expressing GXn in terms of GAn, we obtain: 

Iji := a^E[Gj, trjv(G)] + % - zE[Gji] + E[(GAAr),,] = 0. (4.12) 

Now, we consider the sum ^ U*f.UkjIji, k = 1, . . .r fixed and we denote ak = j Ul^UkjGji = 
{UGU*)kk- Then, we have the following equality, using that U is unitary: 

a^E[ak tYN{G)] + 1 - zE[ak] + ^ ?74C/fc,E[(GAw),,] = 0. 

hi 



Now, 



Y,U*kUkjn{GAN)ji] = E[{UGAMU*)kk]=mUGU*AUU*)kk] 
= 7feE[([/Gf/*)fcfc]=7fcE[afc]. 

Therefore, 

a^E[ak trjv(G)] + 1 + (7fe - z)E[ak] = 0. 
Since ak is bounded and Var(trjv((j)) = 0{^), we obtain 

E[afc](aV(^) +7k-z) + l = O(^). 
Then using (4.10) we deduce that E[ak]{(T'^ga{z) + 7fc - 2:) + 1 = Oijr) and using (3.7) 



Ik 



fe=i 



(4.13) 



(4.14) 



b) The general Wigner case 
We shall prove that (4.13) still holds. We now rely on Lemma 4.2 to obtain the analogue of (4.12) 



Jij := a''E[Gji trjv(G)] + (5^ - zE[Gji] + E[{GAN)ji] + 



K4 

6N 



1 ^ 



1=1 



= 0(^). (4.15) 



The term Aijj is a fixed hnear combination of the third derivative of $ := Gji with respect to Re{XM)u 
(i.e. in the direction en = En + En) and $Jm(X;v)ii (i-e. in the direction fu := yJ^-i{Eii — En)). We 
don't need to write the exact form of this term since we just want to show that this term will give a 
contribution of order 0{j^) in the equation for hN{z). Let us write the derivative in the direction en: 

E[{GeuGeuGeuG)ji] 
which is the sum of eight terms of the form: 

E[Gjii Gi^i^Gi^i^Gi^i] 
where if i2q+i = i (resp. I), then 12^+2 = I (resp. i), q = 0,l, 2. 
Lemma 4.3. Let 1 < k < r fixed, then 

N ^ N 

'TV. 



(4.16) 



N N 

F{N) := I ^ Ut,Ukj-Y,mi,j,i]\ < C\^rnz\ 



(4.17) 



for a numerical constant C. 



Proof: F{N) is the sum of eight terms corresponding to (4.16). Let us write for example the term 
corresponding to ii = i, 13 = i, 15 = i: 



]Y XI ^ikUkjElGjiGuGuGii] 

]y ^ifc fci Gil 

i,l 

^Y.^UUG)ki{G^G''G^) 



E 



E 



where the superscript T denotes the transpose of the matrix and is the diagonal matrix with 
entries Gu. From the bounds ||(jjv(-2)|| < |Sjm2|~^ and \\U\\ = 1, we get the bound given in the 



lemma. 

We give the maj oration for the term corresponding to ii = I, is = I, = I: 



N 



= E 



Its absolute value is bounded by E ^ - ; \Gu ^ |9m2;| ^ and thanks to lemma 3.2 by |9rn2:| The 
other terms arc treated in the same way. □ 

As in the Gaussian case, we now consider the sum J2i,j ^ik^kjJji- From Lemma 4.3 and the bound 
(using Cauchy-Schwarz inequality) 

N 

J2 \u:m<n, 

we still get (4.13) and thus (4.14). More precisely, we proved 
Proposition 4.3. Let h^iz) = E[Tr(^;vGjv(2;))], then 



\hNiz)-J2 



Ik 



4.1.4 Convergence of E 



^^-C7V(^)-7fc 



(at Sfe=i ^fefe) 



<^i&l£n(i. + N). 



We now study the last term in the master inequality of Theorem 4.1. For the non Deformed Wigner 
matrices, it is shown in [K-K-P] that 



Rn{z) :=E 



N >(x> 



Moreover, Proposition 3.2 in [C-D], in the more general setting of several independent Wigner matrices, 

gives an estimate of |i?jv(z) — g^(z)\. The above convergence holds true in the Deformed case. We 
just give some hints of the proof of the estimate of \Rn{z) — g'^{z) \ since the computations are almost 
the same as in the non Deformed case. Let us set 



1 ^ 



We start from the resolvent identity 

zGkk = 



N 



^Ik 



1 + Y,{MN)kiGi, 
1=1 

N N 

1 + ^(^iv)fciGifc + y^(^jv)fc;G;fc 



and 



1 ^ 1 ^ 1 ^ 

zdN{z) = — ^ Gkk + y^X^N^)kkGkk + jy X/ i^N)klGlkGkk- 

k=l k=l k,l=l 

For the last term, we apply an integration by part formula (Lemma 4.2) to obtain (see [K-K-P], [C-D]) 



E 



k,l=l 



{XN)kiGikGkk = cr^E (— Gkk)dN{z) 



fe=l 



It remains to see that the additional term due to is of order 0(-^). 



N r 
1 V— ^ 1 , , 

fe=i p=i 



and 



\^J2(^NG)kkGkk\<{J2\^p\) 
fe=l p=l 

We thus obtain (again with the help of a variance estimate) 



zE[dN{z)] = gN{z) + a^gN{z)E[dN{z)] + 0{^). 
Then using (4.10) and since d^iz) is bounded we deduce that 

1 , 



zE[dNiz)] = ga{z) + cj^g„{z)E[dN{z)] + 0{-^, 



Thus (using (3.7)) 



E[d^{z)]= ^''['\ , +0{h . =gl{z). 

z-cr^ga{z) ^N^N — yoo z - a^gaiz) " 

Now, using some variance estimate, 

E[dl{z)] = m^{z)]f + 0(1) = g'„{z) + 0(1). □ 

We can now give our final master inequality for gN{z) following our previous estimates: 
Theorem 4.2. For z G C such that 9to(z) > 0, gwiz) satisfies 



(T^g%{z) - zgN{z) + 1 + ^E„{z) 



where E^iz) = z-a^g'\z)--yk ^dai^)' fourth cumulant of the distribution (x. 

Note that E„{z) can be written in terms of the distinct eigenvalues 9j's oi An as 

J 



Let us set 

La{z) = gM)~Miz - sr^)E^{z) (4.18) 
where s is a centered semicircular random variable with variance ct^. 

4.2 Estimation of \ga{z) — gN{z) + j;jL„{z)\ 

The method is roughly the same as the one described in Section 3.6 in [C-D]. Nevertheless we choose 
to develop it here for the reader convenience. We have for any z in O^, by using (4.9), 



\ga{z)-gN{z) + jjLaiz)\ 



= |E 
< IE 



{z - s)-\Ki,{z) - s)-\Kn{z) -z) + -ga{z)-\z - s)-^E,{z) 
{z - s)-\Kj,{z) - s)-\Km{z) - z + ^g,{z)-'E,{z)) 



+E [\{z - s)-^{{z - s)-^ - (A^iz) - s)-'}\] -\g,{z)-^EAz)\ 

< 2\^m{z)\~^\AN{z) - z + lE„{z)g„{z)-'\ 

+ \^n{z) - z\{\z\ + K) 

where we made use of the estimates (3.5), (4.8), Vz G C\M, x € M, 

1 



< \'im{z)\-\ 



z — X 

\Ea{z)\ < Pii\^miz)\-') (using (3.7)). 

Let us write 

\An{z) - z + jfE„{z)g„{z)-'\ 

= JiT) W9%{z) - zgNiz) + 1 + jfE^iz)} + ^-^^^±{g^{z) - g,{z)} E,{z). 
We get from Theorem 4.2, (4.6), (4.10), (4.19), (3.5) 

|A^(.) - . + ]^E^{z)g.{z)-'\ < {\z\ + ;^)3 J^i5(iy^) _ 
Finally, using also (4.7), we get for any z in O^, 



Now, for z ^ such that Qm{z) > 

Pe{\{Qm{z)r'\) 



1 < 4- 



< i\z\+K) 



N 

P8(|9m(z)-i|) 



N2 



i\z\+a^\^m{z)\-^)\^m{z)\-^ 



N 



We get 



\gaiz) - gN{z) + ^L„iz)\ < \g„{z)-gN{z)\ + ^\L„iz)\ 

< (|,|^^)W.)n^ 



{\z\+K) 



N 

Pg{\^m{z)\-') , 1 



(4.19) 



N 



+ -P,{\Qm{z)r){\z\+K) 



Thus, for any z such that SJto(z) > 0, 



2Pi7(|9m(z)|-^) 



7V2 



\g4z) - g.iz) + li..(.)| < (N + i,)3^I&|£)n . 



(4.20) 



Let us denote for a while gN = gf/^ and = L^" ■ Note that we get exactly the same estimation 
(4.20) dealing with —An instead of An- Hence since gaiz) = —ga{—z), ^ {^) = ~9n^{~^) (using 
the symmetry assumption on /i) and L~"^"(z) = L^^{—z), it readily follows that (4.20) is also valid 

for any z such that 9mz < 0. In conclusion, 

Proposition 4.4. For any z G C\M, 

\9Az) - 9n{z) + ^L^{z)\ < i\z\ + KfM^^^r^. (4.21) 
4.3 Inclusion of the spectrum of Mjv 

The following step now consists in deducing Proposition 4.6 from Proposition 4.4 (from which we 
will easily deduce the appropriate inclusion of the spectrum of Mat). Since this transition is based 
on the inverse Stieltjes transform, we start with establishing the fundamental Proposition 4.5 below 
concerning the nature of L^. Note that one can rewrite as 



L,{z) = g^z)-' X gUz) [ J^fc, + ) • (4-22) 



We recall that J+a (resp. J-a) denotes the number of j's such that 6j > a (resp. 9j < —a). As in 
the introduction, we define 

2 

pe, =(^3 + Y, (4-23) 
which is > 2(7 (resp. < — 2cr) when Oj > a (resp. < —a). 

Proposition 4.5. L„ is the Stieltjes transform of a distribution A^. with compact support 

Ka{6i,--- := [pejipej.r-,--- ; Pej-j_^+i] ^ [-2a;2a][J ^^pej_^^; pej_^^_r; ■ ■ ■ 

As in [S] , we will use the following characterization: 
Theorem 4.3. [T] 

• Let A be a distribution on M. with compact support. Define the Stieltjes transform of A, I : 
C\R -^C by 

l{z) = a(— 

\z — X ^ 

Then I is analytic in C\M and has an analytic continuation to C\supp{A). Moreover 
(ci) l{z) as \z\ oo, 

(C2y) there exists a constant C > 0, on n G N and a compact set K cM. containing supp{A) such 



that for any z G 

\l{z)\ < Cmax{dist{z,K)-'',l}, 

(03) for any (j) G C°°(M, R) with compact support 

A{(f>) = lim 9m / 4>{x)l{x + iy)dx. 

• Conversely, if K is a compact subset o/M and if I : <C\K C is an analytic function satisfying 
(ci) and (ci) above, then I is the Stieltjes transform of a compactly supported distribution A on 
K. Moreover, supp{A) is exactly the set of singular points of I in K. 

In our proof of Proposition 4.5, we will refer to the following lemma which gives several properties on 



Lemma 4.4. g^- is analytic and invertible on C\[— 2ct, 2ct] and its inverse satisfied 

Za{9) = ^ + <7^5, V5 e 5^(C\[-2a,2a]). 
(a) The complement of the support of is characterized as follows 

X e ]R\[-2(T, 2c7] <s=^ 3g G R* such that \-\ > a and x = Za{g)- 

(h) Given x e M\[— 2(t, 2(t] and G R such that \9\ > a, one has 



This lemma can be easily proved using for example the explicit expression of (derived from (3.3)) 
namely for all x G M\[— 2cr, 2a], 



Proof of Proposition 4.5: Using (4.22), one readily sees that the singular points of L„ is the set 

[— 2(t; 2cr] U 1^ € 2cr, 2cr] and € Spect(Ajv)|. Hence (using point (b) of Lemma 4.4) the set 
of singular points of L„ is exactly K„{di, ■ ■ ■ , 6j). 

Now, we are going to show that L^- satisfies (ci) and (02) of Theorem 4.3. We have obviously that 

\z-a^g,{z)-e,\>\ \z-ej\-\a^g„{z)\ \. 

Now, let a > such that a > 2cr and for any j = 1, . . . , J, a — > For any 2; G C such that 

\z\ > a, 

\z-eA>\z\-\0A>-^ 

a — za 

and according to (3.8) 



K9a{z)\ < , , , < 



\z\ —2a a — 2a 

Thus we get that for 2; e C such that \z\ > a, 

2 

\z-a'g^iz)-e,\>\z\-\9j\-^^. 
Using also (3.8), (3.9), (3.10), we get readily that for \z\ > a, 

Then, it is clear than |i(^(2;)| — » when \z\ +00 and (ci) is satisfied. 

Now we follow the approach of [S] (Lemma 5.5) to prove (02). Denote by £ the convex envelope of 
Ka{9i,- ■ ■ ,0j) and define the interval 

K ■.= {xG R;dist{x,£) < 1} = [min{a; e K^iei,--- ,6»j)} - l;max{a; e K^iOi,--- + l] 

and 

D = {zeC;0< dist{z, K) < 1} . 



• Let z e D n C\M with ^e{z) G K. We have dist{z,K) = |3mz| < 1. Using the upper bounds 
(3.4), (3.5), (3.6) and (3.7), we easily deduce that there exists some constant Co such that for 
any 2 e £1 n C\M with ^e{z) € K. 

\La{z)\ < Co\'^mz\-'^ = Codist{z,K)-'^ = Co max(rfisi(2, Jf)"^; 1). 

• Let z e Dn C\R with 3fJe(z) ^ K. Then dist{z, K^{0i, ■■■ , Oj)) > 1. Since is bounded on 
compact subsets of 'C\K„{9i, • • • , 9j), we easily deduce that there exists some constant C\ such 
that for any z & D with ^e{z) ^ K, 

\La{z)\ < Ci < Cidist{z,K)-'^ = Cima,x{dist{z,K)-'^; 1). 

• Since |i<7(^)| ~* when \z\ — > +oo, is bounded on C\D. Thus, there exists some constant 
C2 such that for any z G C\D, 

\La{z)\ < C2=C2max{dist{z,K)-'^;l). 

Hence (02) is satisfied with C = max(Co, C\,C2) and n = 7 and Proposition 4.5 follows from Theorem 
4.3. □ 

We are now in position to deduce the following proposition from the estimate (4.21). 
Proposition 4.6. For any smooth function (p with compact support 

E[tvMi^{MN))] = y^d/x,e + ^Aa(v)+0(^). (4.24) 

Consequently, for (p smooth, constant outside a compact set and such that supp{ip)C\K^{6i, • • • , 9j) = 

9, 

tiNi^iMN)) = 0{^) a.s. (4.25) 

iV 3 

Proof: Using the inverse Stieltjes tranform, we get respectively that, for any (p in C°°(M, M) with 
compact support, 

E[trjv(</?(Mjv))] - / ipdiJ.sc - ^^^^ = -- lim 9m / ip{x)rN{x + iy)dx (4.26) 

J JM IT y^0+ 7r 

where rjv = 9a{z) — 9n{z) + jfL„{z) satisfies, according to Proposition 4.4, for any z G 



\rN{z)\ < -^{\z\+KrPk{\'^m{z)-'\) (4.27) 

where a = 3 and k = 17. 

We refer the reader to the Appendix of [C-D] where it is proved using the ideas of [H-T] that 

lim sup I / ip{x)h{x + iy)dx\ < C < +00 (4.28) 

y^o+ Jm. 

when h is an analytic function on C\R which satisfies 

\h{z)\ < {\z\ + KrPk{\^m{z)-'\). (4.29) 
Dealing with h{z) = N'^rj\f{z), we deduce that 



f C 
limsupl / ip{x)rN{x + iy)dx\ < (4.30) 



*0+ 



and then (4.24). 

Following the proof of Lemma 5.6 in [S], one can show that Acr(l) = 0. Then, the rest of the proof of 



(4.25) sticks to the proof of Lemma 6.3 in [H-T] (using Lemma 3.1). □ 

Following [H-T](Theorem 6.4), we set K = KAOi.- ■ ■ , 6',7) + (-|, §), F = {teW.; d{t, K„{9i,- ■ ■ , Oj)) > 
e} and take (p e C°^(M,K) such that < <y9 < 1, ip{t) = for t e if and ip{t) ^ 1 ioT t e F. Then 
according to (4.25), trjv(<^(Mjv)) = O(^) a.s. Since > If, it follows that trjv(lF(Afjv)) = 0{^) 

N'S N'S 

a.s. and thus the number of eigenvalues of Mjv in F is almost surely a O(^t-) as N goes to infinity. 

Since for each N this number has to be an integer we dediice that the number of eigenvalues of Mn 
in F is zero almost surely as N goes to infinity. The fundamental inclusion (2.16) follows, namely, for 
any e > 0, almost surely 

Spect(MN) C K^(6li,--- ,6ij) + (-£,£) 

when N goes to infinity. 

Such a method can be carried out in the case of Wigner real symmetric matrices; then the approx- 
imate Master equation is the following 

a^g^[zf-zg^[z)+l+-!^E[{- J2 Gkkizrf] (tr^ Giv(^)')+E (tr^ [A^G^(z)]) = O(^). 

fc=i 

Note that the additionnal term (tr^r Gn{z)^) already appears in the non-Deformed GOE case in 
[S]. One can establish in a similar way the analog of (4.10) and then, following the proof of Corollary 
3.3 in [S], deduce that 

E (trjv GN{zf) =M{{z- s)-2) + 0(1), 

where s is a centered semi-circular variable with variance a^. Hence by similar arguments as in the 
complex case, one get the master equation 

a^QNizf - zg^iz) + 1 + ^E,{z) = O(^) 

where 

It can be proved that Lcr{z) := grr{z)~^'E{{z — s)~'^)E„{z) is the Stieltjes transform of a distribution Acr 
with compact support K„{Oi,- ■ ■ , 9j) too. The last arguments hold likewise in the real symmetric case. 

Hence we have established 

Theorem 4.4. Let J+o- (resp. J-„) he the number ofj's such that 6j > a (resp. Oj < —cr). Then for 
any £ > 0, almost surely, there is no eigenvalue of Mjv in 

(-0O, P0J - e) U {p0j + e, p0j_, - e) U • • • U (pe^. + e, -2a - e) 
U {2a + e,p0j^^-e)[J---[J{p0^+e,p0,-e)[J{p0,+e,+oo) (4.31) 

when N is large enough. 

Remark 4.1. As soon as e > is small enough that is when 

2e < min (pe^^i - Pe,,J - J-a + 2<j<J or 2 < j < J+^; -2a - p0j_j_^^^ ; p0j^^ - 2a) 
the union (4-31) is made of non-empty disjoint intervals. 



4.4 Almost sure convergence of the first extremal eigenvalues 

As announced in the introduction, Theorem 2.1 is the analog of the main statement of [Bk-Sl] es- 
tablished for general spiked population models (1.1). The previous Theorem 4.4 is the main step of 
the proof since now, we can quite easily adapt the arguments needed for the conclusion of [Bk-Sl] 
viewing the Deformed Wigner model (1.2) as the additive analog of the spiked population model (1.1). 

Let us consider one of the positive eigenvalue 9j of the Aj^^a. We recall that this implies that 
Afei+...+/c^_j+i(Ajv) = 9j for all 1 < z < kj. We want to show that if 6*^ > a (i.e. with our notations, if 
j G {1, • • • , J+cr}), the corresponding eigenvalues of Mjv almost surely jump above the right endpoint 
2a of the semicircle support as 

yi<i<kj, Afei+...+fe^._i+i(Mjv) — > pOj a.s. 

whereas the rest of the asymptotic spectrum of Mjv lies below 2a with 

>'ki+-+kj^^+i{MN) — > 2a a.s.. 

Analog results hold for the negative eigenvalues 9j (see points (c) and (d) of Theorem 2.1). To describe 
the phenomenon, one can say that, when A'' is large enough, the (first extremal) eigenvalues of Mn 
can be viewed as a "smoothed" deformation of the (first extremal) eigenvalues of An- So, our main 
purpose now is to establish the link between the spectra of the matrices Mjv and A^. According to 
the analysis made in the previous section (Proposition 4.5), we yet know that the 6j's are related to 
the poj's through the Stieltjes transform g^. More precisely, one has 

1 

for all j such that \6j\ > a, 



Actually, one can refine this analysis and state the following important Lemma 4.5 on As before, 
we denote (recall Lemma 4.4) by its inverse which is given by 

Using Lemma 4.4, one readily sees that the set '^K„{6i,- ■ ■ , 6j) can be characterized as follows 

X e "K^iOi,- ■■ Jj) <^ 3g e such that x = z„{g) (4.32) 

where 



:= jgr e M* : |-| > fT and - ^ Spect(Aiv)}. 
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Obviously, when x G '^Ka{9i^ ■ ■ ■ ,0j), one has g = ga{x). 

Lemma 4.5. Let [a,b] he a coiv,pact set contained in '^Kc{9\, - ■ ■ ,9j). Then, 

9a{a) 9a{b) 

(ii) For all < a < a, the interval [za{9cr{o))j z^{ga{b))] is contained in '^K^{6i,--- ,6j) and 
za{9a{b)) - za{ga{a)) >b-a. 

Proof: The function being increasing, (i) readily follows from (4.32). 

Noticing that Qc C Ga for all ct < a implies (recall also that ^o- decreases on [a, b] ) that [ffo- (6) ; 5^ (a)] C 
Ga- Relation (4.32) combined with the fact that the function za is decreasing on [gaib); gaia)] leads 
to 

[zaigAa));zaigAb))] C ^Ka{9i,- ■ ■ ,0j) 
and the first part of (ii) is stated. Now, we have 

laia) := za{ga{b)) - za{ga{a)) 

^ ^ +er''{g„{h)-g„{a)). 



9a{b) 9Aa) 



Since decreases on [a, 6], we have ga{i>) — ga{<i) ^ and thus la is decreasing on M+. Then the last 
point of (ii) follows since = b — a. □ 

Thanks to this lemma and the previous Theorem 4.4, one can state the asymptotic relation between 
the spectrum of Ajv and the one of Mj^. 

Let [a, b] be an interval contained in ''K„{6i, ■ ■ ■ , 9j). By Theorem 4.4, [a, 6] is outside the spectrum 
of Mjv. Moreover, from Lemma 4.5 (i), it corresponds an interval [a' , b'] outside the spectrum of Aj^r 
i.e. there is an integer in G {0, . . . , N} such that 

Ai„+i(^jv) < := a' and Xi^iA^) > := b' . (4.33) 

a and a' (resp. b and b') are linked as follows 

a = Pa' ■■= a + — (resp. b = pb'). 

Our aim now is to prove that [a, b] splits the eigenvalues of Mjv exactly as [a', b'] splits the spectrum 
of ^jv- In [B-S2], one talks about the exact separation phenomenon. 

Theorem 4.5. With iN satisfying (4.33), one has 

P[Ai^+i(Mjv) < a and Ai„(Mjv) > b, for all large N] = 1. (4.34) 

Remark 4.2. This result is the analog of the main statement of [B-S2J (cf. Theorem 1.2 of [B-S2]) 
established in the spiked population setting (and in fact for quite general sample covariance matrices). 

Intuitively, the statement of Theorem 4.5 seems rather natural when a is close to zero. Indeed, 
when N goes to oo, since the spectrum of -^Wn is concentrated in [—2a, 2a] (recall (2.9)), the 
spectrum of would be close to the one of A^ as soon as a will be close to zero (in other words, 
the spectrum of Mjv is, viewed as a deformation of the one of An, continuous in ct in a neighborhood 
of zero). Actually, this can be justified regardless of the size of a thanks to the following classical 
result (due to Weyl). 

Lemma 4.6. (cf. Theorem 4-3.7 of [H-J]) Let B and C be two N x N Hermitian matrices. For any 
pair of integers j, k such that 1 < j,k < N and j + k < N + 1, we have 

Xj+k-i{B + C) < Xj{B)+Xk{C). 

For any pair of integers j, k such that I < j,k < N and j + k > N + 1, we have 

Xj{B) + Xk{C) < Xj+k-N{B + C). 

Note that this lemma is the additive analogue of Lemma 1.1 of [B-S2] needed for the investigation of 
the spiked population model. 

Proof of Theorem 4.5: With our choice of [a, b] and the very definition of the spectrum of the 
j4jv's, one can consider e' > small enough such that, for all large N, 

Ai„+i(Ajv) < -^r^ - e' and Ai„(Ajv) > -^jtt + e'- 
9cr{a) ga{b) 

Given L > and fc > (their size will be determined later), we introduce the matrix W'^^ = 
-Wn and M^'^ = An + -j^W^^ . We also define 



c7fe,L = , cr, ak,L = Adcria)) and bk,L = z^^ L{9a{b)) 



where we recall that Zq-j. ^(5) = ^ + ffe i.9- Note that for all L > 0, one has: = Mn, ao,L = cl 

and 60, L = b. 

We first choose the size of L as follows. We take Lq large enough such that for all L > Lq, 

max(la2(|54a)| + |5.WI);^) <^ (4-35) 

From the very definition of the a^^i's and 6fe,L's, one can easily see that bk,L — <ik,L > b — a (using 
the last point of (ii) in Lemma 4.5) and that this choice of Lq ensures that, for all L> Lq and for all 
fc > 0, 

I , b — a ,1, , I b — a , , 

|afc+i,L - afc,L| < and \bh+\,L - Ofc.il < (4.36) 

Now, we fix L such that L > Lq and we write = ak,L, bk = bk,L and cr/c = ak,L- 

Lemma 4.6 first gives that 

Xi^+i{M^'')<ak-e' -alg^{a) + ^^=Xi{^WN) for In < N 

and X,^{M^^)>bk + e' -alg^{b) + ^^=XN{^WN) for iN > 0. 

Furthermore, according to (2.9), the two first extremal eigenvalues of -^Wn are such that, almost 
surely, at least for N large enough 

< max(-Ajv(^VFAr), Xi{^Wn)) < 3a. 

y N y N 

Thus, for all fc, almost surely, at least for N large enough (A'' does not depend on k), 

< J- X max(-AAr(^VFAr),Ai(^W^jv)) < 3afe. 

1 -I- 1. V iV V-/V 

As C7fe — > when k +00, there is K large enough such that for all fc > if, 

max(|3<Tfe - algcr{a)\, |3c7fe + (jlga{b)\) < e 
and then a.s for N large enough 

Ai„+i(M^'^) < ak if IN < N, (4.37) 

Ai„(M^'^) >6fc ii iN>0 (4.38) 

Since (4.37) respectively (4.38) are obviously satisfied too for ij^ = N resp. Zn = 0, we have established 
that for any in & {0, . . . , N}, for all k > K, 

P[Ai„+i(M^'^) < ak and Ai^(M^'^) > bk for aU large N] = 1. 

In particular, 

F[Xi^+i{M^'^) < GK and Xi^{M^'^) > for all large N] = 1. (4.39) 

Now, we shall show that with probability 1: fov N large, [aK,bK] and [a,b] split the eigenvalues 
of, respectively, M^'^ and Mjv having equal amount of eigenvalues to the left sides of the intervals. 
To this aim, we will proceed by induction on k and show that, for all fc > 0, [ak,bk] and [a,b] split 
the eigenvalues of M^^ and Mjv (recall that Mn = M^^) in exactly the same way. To begin, let us 
consider for all fc > 0, the set 

Efc = {no eigenvalues of M^^ in [uk, bk], for all large A''}. 



By Lemma 4.5 (ii) and Theorem 4.4, we know that P(Efc) = 1 for aU k. In particular, from the fact 
that P(Eo) = 1, one has for all ui € Eq, for all large N: 

3jN{uj) e {0, • • • , A/"} such that Aj„(^)+i(Mjv) < a and Aj„(^)(Mjv) > b. (4.40) 

Extending the random variable Jn by setting for instance Jn := — 1 on '^Eq, we want to show that 
for all k, 

P[A,v+i(M^'^) < ak and A,v(M^) > bk, for all large N] = 1. (4.41) 

This can be done by induction calling, once more time, on Lemma 4.6. By (4.40), this is true for 
fc = 0. Now, let us assume that (4.41) holds true. We shall show that this still holds replacing k by 
k + 1. One has 

fk+hL _ Mk,L .1 1.1 



A,-,+i(M^+i'^) < A,,+i(Mr ) + 1(-A^(-^T^^)). 



so, by Lemma 4.6, 



But, for A'' large enough, < —Xn{-^Wn) < Sa a.s, so by the choice (4.35) on L, 



Similarly, one can show that 



A,-^+i(M^+^'^) <ak+ ^—^ := afe. 



a.s. Aj„(M^+^'-^) >hk- ^—r^ := h- 



Now, by (4.36), one readily observes o/j — Uk+i = ak — afc+i + > and similarly that bk — bk+i < 
which implies that 

[ak,bk] C [afc+i,5fc+i]. 
So, as P(Efe_|_i) = 1, we deduce that with probability 1: 

Aj„+i(M^+^'^) < Qk+i and \j^,{M^'^^'^) > bk+i, for all TV large enough. 

As a consequence, (4.41) holds for all A; > and in particular for k = K. Comparing this with (4.39), 
we deduce that Jn = in a.s. and 

P[Ai„+i(MAr) < a and Ai„(MAr) > b for all large N] = 1. 

This ends the proof of Theorem 4.5. □ 

Now, we are in position to prove the main Theorem 2.1. 

Proof of Theorem 2.1: Our reasoning is close to the last Section 4 of [Bk-Sl]. It is enough to 
establish parts (a) and (6) since the assertions (c) and {d) can then be deduced by taking -Mn instead 
of Mjv. 

The proof of (a) is mainly based on successive applications of Theorem 4.5. Fix an integer 1 < j < J+a, 
and let us consider for e > 0, the interval [a,b] = [pg^ + e,pg._-^ — e] which is included in the union 

(4.31) (at least for e small enough). We define Kj = ki-\ h kj (with 9o +oo and the convention 

Xo{Mn) = \o{An) = +oo). Since l/g„{pQ^) = 6k ioi k = j — 1 and j and since the function 1/ga is 
continuous and increasing on [a,b], the compact interval [aj)] satisfies (4.33) with ijv = Kj-i (with 
the convention that ijv = if j = 1). Hence, by Theorem 4.5, one has 

P[Aif,_i(Mjv) > peJ_^ - e and AK,_i+i(MAr) < pe. + e, for N large] = 1. 

Similar arguments imply that, for all j e {1, ... , J+a ~ 1}) 



P[Ak, (Mjv) > pej - e and Ajf,+i(Mjv) < pej+, + e, for N large] = 1. 



As a result, we deduce that for all 1 < j < J+a — Ij 

PK - e < Ak, {Mn) <■■■< Ak,_i+i(Mjv) < pej + e, for N large] = 1. (4.42) 

So, letting e go to zero, we obtain (a) for each integer j of {1, • • • , J+o- — !}■ 

Let us now quickly consider the case where j = J+a- Note first that, from the preceding discussion, 
we still have (for e small enough) 

¥[\kj^^_,+^{Mm) < P0j^^ + e, for N large] = 1. 

Then, using the fact that I/Qo- increases continuously on ]2a, +oo[ with l/g„(\2a, +oo[) =]a, +oo[, one 
can show that once e > is small enough, the compact set [a, b] = [2a + e, POj^^ — e] satisfies the 
assumptions of Theorem 4.5 with iN = Kj^^ . This leads to 

^Xkj^^ (Mn) > P9j^^ - € and Xkj^^+i{Mn) < 2c7 + e, for N large] = 1. 

Thus, letting e ^ 0, we deduce that (4.42) holds for j = J_^_„ and the assertion (a) is established. For 
point (6), it remains to prove that 

liininf Aifj^_^+i(MAf) > 2a a.s. 

Such an inequality follows from the fact that the spectral measure of converges almost surely 
towards the semicircle law which is compactly supported in [—2a, 2a]. This completes the proof 
of Theorem 2.1. □ 

5 Fluctuations 

The (complex or real) Wigner matricial models under consideration are the same as previously (i.e 
defined by (i) in Section 2) but now we assume that the perturbation An is diagonal with unique 
eigenvalue > a: = diag(0, 0, • • • ,0). According to the previous section, the a.s convergence of 
Ai(Mjv) towards pg is universal in the sense that it does not depend on /j,. In this section, we are 
going to show that the fluctuations of Xi{Mn) around this universal limit are not universal any more. 
Indeed, we are going to prove that •\/iV(l — fj-)~^(Ai(M]v) — pg) converges in distribution towards 
the convolution of p, and a Gaussian distribution. Hence, the limiting distribution clearly varies with 
p and in particular cannot be Gaussian unless p is Gaussian. 

5.1 Basic tools 

We start with the following results which will be of basic use later on. Note that in the following, a 
complex random variable x will be said standardized if E(a;) = and E(|a;|^) = 1. 

Theorem 5.1. (Lemma 2.7 [B-Sl]) Let B = {bij) be a NxN Hermitian matrix andY^ he a vector of 
size N which contains i.i.d standardized entries with bounded fourth moment. Then there is a constant 
K > Q such that 

EIY^BYat - TrB|2 < KTv{BB*). 

Theorem 5.2. (cf. [B-Ya2] or Appendix by J. Baik and J. Silverstein) Let B = {bij) be a N x N 
random Heiinitian matrix and Yn = (jji, ■ ■ ■ , Vn) be an independent vector of size N which contains 
i.i.d standardized entries with bounded fourth moment and such that E(j/^) = i/ j/i is complex. 
Assume that 

(i) there exists a constant a > (not depending on N) such that \ \B\ \ < a, 
(ii) ^TrB2 

converges in probability to a number 02, 
r*^*/' TJ Sj^i ^li converges in probability to a number a^. 



Then the random variable -^{Y^BYn — TiB) converges in distribution to a Gaussian variable with 
mean zero and variance 

iE\yi\*-l-t/2)ai + {t/2)a2 
where t = A when Yi is real and is 2 when yi is complex. 

Proof: This result is in fact a particular case of a more general result of [B-Ya2] (Theorems 7.1 and 
7.2) which follows from the method of moments. We give an alternative elegant proof by J. Baik and 
J. Silverstein in the Appendix of the present paper. 

Theorem 5.3. (Theorem 1.1 in [B-Yal]) Let f be an analytic function on an open set of the complex 
plane including [—2a;2a]. If the entries {{WN)ij)i<i<j<N of a general Wigner matrix Wn satisfy the 
conditions 



n\{WN)ij\^) = M forii^j, 

For any rj > 0, limAr^+<x, Eij ^ l(^Jv)y |''l{|(w^,)y|>r,V^ 



= 0, 



then N{tYN{fi^^WN)) — j fd^Xgc} converges in distribution towards a Gaussian variable. 

In our setting, jj, satisfies a Poincare inequality and thus, as already noticed in Section 2, fi satisfies 
J \x\'^dfi{x) < +00 for any q in N. Hence, the general Wigner matrices we consider obviously satisfy 
the conditions of Theorem 5.3. Nevertheless, in the following study of fiuctuations, we do not use 
the Poincare inequality; thus one can expect that Theorem 2.2 is still valid under assumptions on the 
only first moments of provided one can prove the a.s convergence of Xi{Mn) towards pe under these 
weaker assumptions. 

5.2 Proof of Theorem 2.2 

The approach is the same for the complex and real settings and is close to the one of [P] and the ideas 
of [B-B-P]. Let Mm-1 be the N — 1 x N — 1 matrix obtained from Mjv removing the first row and 
the first column. Thus, -^S^Myv-i is a non-Deformed Wigner matrix associated with the measure 

II. We will denote by Ai(Mjv^i) the largest eigenvalue and by Ajv-i(Mjv-i) the lowest eigenvalue of 
the matrix Mn-i- Let < S < Let us define the event 



{Ai(iQjv-i) <2a + d; Ajv-i(Mjv-i) > -2ct - <5; Ai(Mjv) > - <5} . 



According to [B-Yi] and Theorem 2.1 of the previous section, limAr^+co IP(f^Af) = 1- Thus, it is 
sufficient to restrict ourselves to the event flN in order to study the convergence in distribution of 
\/]V(Ai(M;v)-pe). 



Let V = ^ {vi, . . . , vn) be an eigenvector corresponding to Xi{Mm)- 
MnV = Xi{Mn)V 



{MN)iiVi+J2f=2iMN-i)ijVj = Xi{MN)v^, yi = 2,...,N. 
Define the following vectors in C^~^: 

V= *(%,..., ujv) 

and 

Ml = * {{Mm)2i, {Mm)ni) = * ((VF^)2i, ■ ■ ■ , {Wn)ni) ■ 

On Qn, Xi{Mn) is not an eigenvalue of Mn-i and one can write the eigen-equations using the 
resolvent G{Xi{Mn)) := {Xi{Mn)In-i - Mjv-i)"^ as follows: 

V = viG{Xi{Mm))M.i. (5.1) 

Xi{Mn)vi = evi + VI + viM\G{Xi{Mn))M.i. (5.2) 



Since vi is obviously non equal to zero, one gets from (5.2) 

Xi{Mn) = e+ + M.*iG(Ai(Miv))M.i. 



(5.3) 



On fijv, Pe is not an eigenvalue of Mjv-i and the resolvent G{pe) := (p^/jv-i — Afjv-i) ^ is well 
defined too. Thus, (5.3) is equivalent to 

M{Mn) ~Pg= + M\G{pe)M.i - ^ + M* \g{\i{Mn)) - G(/Je)l Mi- (5.4) 



Using G{Xi{Mn)) - G{p0) = -(Ai(Mjv) - pe)G{pe)G{Xi{MN)) and Sf<,(p0) = i, one gets 
Xi{MN)-pe 



(Wn) 
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/N 
+M.1 

Finally, defining /6/(z) 
ing identity on Ojv: 



M*,G{pg)M.i-a^gM 



-(Ai(Mjv) - Pe)G{pe) {G{pe) - {Xi{Mn) - pe)G{pe)G{Xi{MN))) 



Ml. 



pe- z 



l|z|<2<T+<5, we can easily deduce from the previous equality the follow- 



{1 + CN + di{N) + d2{N)} VN{Xi{Mn) - pe) = (W^jv)ii + 



N 



N-1 



N 



N-1 



where 



CN = (T trN-1 



fe{MN-i) 



dN = VN^{m*, (g(p^)1||S?„_^||<2, J Ml - ahrM-i (g(p,)] 



\MN-i\\<2a+5 



)} 



5i(iV) = -(Ai(Mw) - pe)M*i GM G(Ai(MAr))M.ila, 



52(iV)=M.* G{p0)l. 



Mjv-i||<2(7+i 



Ml - a'^trN-1 



G{pe)l, 



Mjv-i||<2(7+6 



First 



SsiN) = a^y/N^l^tTN-i (/e(M;v-i)) - J fedfXsJj 

MN)\ < |Ai(Mjv)-pe|||Mi||2||G(pe)||'||G(Ai(M;v))||la„ 

1 



1 ^ 1 

i=2 J 



(P0 - 2(7 - 2S){p0 - 2(7 - 5)2 I A/- ^_ 

(using point (w) of Lemma 3.2). By the law of large numbers jfJ2f=2 converges a.s. to- 

wards (7^ and according to Theorem 2.1, |Ai(Mjv) — pe\ converges a.s. to zero. Hence Si{N) converges 
obviously in probability towards zero. 

Now, since fe is analytic on an open set including [—2(7; 2(7], we deduce from Theorem 5.3 the con- 

/2 
fgdp^sc = 02 2 • 

According to Theorem 5.1 and using Lemma 3.2(v), 

-E t 



trjv 



< 



< 



N-1 \ 

^e(||G(p,)|| 



G{pe)\ 



Mn_i||<2<7+i5 



N 



||Mjv-i||<2<7+5 



K 



TV - 1 (pe - 2(7 - 5Y ■ 



The convergence in probability of 62{N) towards zero readily follows by Tchcbychev inequality. 
Let us check that G{p$)l^^j^^ i\\<2ct+5 satisfies the conditions of Theorem 5.2. 

(i) l|S(P«)l||M„_i||<2a+5ll ^ pe-L-s by Lemma 3.2 (v). 

(ii) According to Theorem 5.3, trjv_i /g (Mjv-i) converges in probability towards / fgdjiac- Since on 
the event {||Mjv-i|| < 2(j + 5}, with limiting probability 1, trjv-i[G(pe)l||j^^ i||<2<7+(5]^ coincide with 
trjv-i /g (Mjv-i), it also converges in probability towards J fgdnsc- 

(iii) It is proved in Proposition 3.1 in [C-D] that for any z e C such that Qmz > 0, X^ila^d'-'l^)]^)^ 
converges in probability towards g'^{z). The same result holds for J2iLi^ i[^{^)]ii)'^^\\M,^ i\\<2<t+6- 
For any e > and any a > 0, 



( 1 ^"^ 

+^ Vjf^l E ([^(^^ + ^")]»)'l||M„_,||<2.+^ ^ 9l{pe + > I j 
+^{\gl{pe)-gl{pe + ia)\>^^. 

Since 

{mP0)U? - {[G{pe + ia)l,f] 1||M„_^,|<2,+, 

= [G{pe) - G{pe + ia)]u[G{pe) + G{pe + ^a)]^^l||M„_l||<2a+^ 
= ia[G(pe)G(pe + ia)]ii[G{p0) + G(pe + ^a)]iil||M„_i||<2a+5' 
we get by using Lemma 3.2 (v) 

\mpe)h? - +*a)].0^|l||M„_.||<2.+. < ^p,-2a-5r 

Similarly, we get that 

\9l{Pe) - glipe + < 



Thus, choosing a such that — _^"_.sii < I, we readily deduce the convergence in probability of 



{pe-2af 

(pe-2<T-5)^ ^ 3' 
N-1 

L 

TV- 



1 ^-1 ^ 

— E(['^(^^)]")^^||MAr-i||<2<7+5 
i=l 



towards g^{p0). 

Since i||<2cr+(5 ^'^^ independent, we can deduce from Theorem 5.2 that rfjv con- 

verges in distribution towards a Gaussian law with mean zero and variance 

4 /I'll?/! (^^)l2 |4n 1 1 * 

ve:=a |(E(|^— I )-l-i/2)^) + -^^ 
where t = 4 in the real setting and f = 2 in the complex one. One readily verifies that ve satisfies (2.12). 



Let < e < 1. Since 6i{N) + d2{N) converges in probability towards zero, the probability of the event 
flj^ = Q]^ n {1^1 (A?') + d2{N)\ < e} tends to 1. Now, since cjv > we have the following identity on 

x/]V(Ai(Mjv) -Pe) = ^ |(W-Ar)ii + \f^d^ + ^ jJ^Ss{N)^ 

with UN '■= I + cn + 6i{N) + 52{N) converging in distribution towards (1 — fr)^^- Moreover, 
since {Wn)ii and cIn are independent, (Wjv)ii + j^zjdN + \J ^^^^^{N) converges in distribution 
towards the convolution of /i and a Gaussian distribution Af{0,vo). Finally, we can conclude that 
y/N (Ai(Mjv) — pe) converges in distribution towards (1 — ^) {/U * A/'(0, vq)}. □ 

Appendix by J. Baik and J. Silverstein 

This Appendix presents the proof by J. Baik and J. Silverstein of the CLT (given by Theorem 5.2) 
needed in the previous section for the proof of Theorem 2.2. Their proof is based on a writing of the 
expression 

{l/y/N){Y^BYN - TrB) (A.l) 
as a sum of martingale difference, and uses the following CLT. 

Theorem A.l. (Theorem 35.12 of [BiJ) For each N, suppose Zni,Zn2, ■ ■ ■ , Z^tn real martingale 
difference sequence with respect to the increasing a -field {Tn,]} having second moments. If asN ^ oo, 

^EiZ%,\J^N,j-i)^v' (A.2) 

where is a positive constant, and for each e > 0, 

J2nZ%jliz.,l>e)^0 (A.3) 

then 

Proof of Theorem 5.2: First, one can write (A.l) as a sum of martingale differences: 

N , 

{1/Vn){Y^BYn - TrB) = (I/Vn) ^ f - 1)6^^ + ViY. Vi^-^ + E 

N , . N 

1=1 ^ j<i j<i ' i=l 

where 

Zi = Zm = (i/\/iv)((i2/i|2 - i)bii +yiJ2 yjbij +yiYl yj^a^- 

j<i j<i 

Let J^jv^i (resp. J^n,o) be the <j-field generated by yi, ... ,yi and B (resp. by B). Let also Ej(-) denote 
conditional expectation with respect to J^N,i- 

It is clear that Zj is measurable with respect to J^N,i and satisfies Ej_i(Zj) = 0. 
We will show the conditions of Theorem A.l are met. 



To verify the Lindeberg condition (A. 3), we need to show this property is closed under addition. 
This will follow from the following fact. For random variables Xi, X2, and positive e 

Ei\X,+X2\^^X,+X.\>e)) < ^E{\X,\H(^\x,\>e/2))+mX2\^M\X,\>e/2))). (A.4) 

Indeed, we have 

Edv'flP i(\Xi+X2\>e)) < l(|Xi|>£/2)) + E(|^lP l(|Xi|<e/2,|X2|>e/2)) 

< E(|Xi|2 l(|x,|>e/2)) + {eVm\X2\ > e/2) 

< E(|Xi| l(|Xi|>e/2)) +E(|X2| l(|X2|>e/2))- 

The same bound starting with X2 leads (A.4). 

Write Zi=Xi+ X^, with = (l/\/]V)(|yi|2 - 1)6^^. Then for e > 0, 

N 

^E(|Xi r \\xii>e)) < a'E{{\y,\^ - if l(||,^p_i|>vive/a)) - (A.5) 

i=l 

as A'' —> 00, by dominated convergence theorem. 
We have 

nT.y^hjf = E(|yi|4 5^|6,,f) + 2E(5^|&,,J>,,.,|2)+E(|y?|2^6?,J?,.J 

j<i j<i * * 

< a'[E\y,\^ + 2 + E\yl\'] 

where the sum ^ is over {ji < i, j2 < i, ji ^ j2}- Therefore EIATII'' = o{N~^) so that for any e > 0, 

* 

N N 

j2n\n\'H\xi\>e))<{y^')J2^\n\'^o bsn^oo. (a.6) 

i=l i=l 

Thus, by (A.5), (A.6) and (A.4), {ZJ satisfies (A.3). 
Now, we shall verify condition (A. 2). We have 

AT N ^ 

^E,_iZf = (l/Ar)^|(E|yi|4-l)62 +Ey1(^y,6,,)2+Ey2(^^,6,,)2 (A.7) 

i=l i=l ^ j<i j<i 

+2E(|yi|2yi)6,, yjb.j + 2E{\yi\'^ yi)hi ^ y^hj + 2(^ %&ij)(X! VjK) 

j<i j<i j<i j<i 

Let Bl denote the strictly lower triangular part of B. We have 

N 

E[{l/N)Y,biiY.yAj]=0 

i=l j<i 

and using Cauchy-Schwarz 

N 



E\{i/N)j2h.J2y^h 



i=l j<i 

N-1 N-1 



= EKl/AT) J2 yjJ2biibij\^ = iW^MYl J2biibijJ2bii^ij) = {l/N^)E{J2biibH{BLBl)ii) 

j=i i>j j=i i>j i>j it 

< E[(max&i,)2(l/n)(^ \{BlBI)u\Y/^] = E[(maxfo,,)2(l/A^)Tr((5i5£)2)i/2] 

ii 

< E[{ma^buf{l/VN)\\BLf]. 



We apply the following bound (due to R. Mathias, sec [Mt]): ||-Bl|| < 7Af||-S||) where = 0(ln7V) 
and the bound ||-B|| < a to conclude that 

N 

as N 00. 

Then (recall that Eyf = when yi is complex), (A. 7) can be written as 

N N 

J2E^-iZf = (l/iV)^[(E|yi|^-l)62.+f(^y.6,.)(^y.6,.)]+o,.p(l) 

i—1 i—1 j<i j<i 

N 

= (l/N) ^(E|2/i|4 - l)bl + t{l/N)Y*BlBLYN + o,.p(l) (A.8) 

where t = A when t/i real, and is 2 when yi is complex. 

Besides, from Lemma 2.7 in [B-Sl] (recalled in Theorem 5.1) we have 

n{l/N){Y;,BlBLYM - TrS£iJi)|2 < (l/iV2)E(Tr(S2Si)2) < Rnsf^^ ^ 
as A'' — > oo. So, as 

TvBlBi^ = ^ |6,,f = (l/2)(TrS2 - 

i<i i 

(A.8) implies that condition (A. 2) holds with 

?;2 = (E|j/i|4-l-t/2)a? + (t/2)a2. 

Thus, by Theorem A.l, we deduce that {1/Vn){Y^BYis[—TiB) converges in distribution to a Gaussian 
variable with mean zero and variance v^. □ 
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